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ABSTRACT 

The existence of black holes (BHs) of mass ^ 10® Mq at z > 6 is a big puzzle in astrophysics because 
even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 Mq to grow 
into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive 
stars with mass M ^ 10® Mq into massive seed BHs. We find that if a jet is launched from the 
accretion disk around the central BH, the jet can break out the star because of the structure of the 
radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^- 
10® s and flux of 10“^^-10“® ergs“^ cm“® could be detectable by Swift. We estimate an event rate of 
< lyr“^. The total explosion energy is > 10®®-10®® erg. The resulting negative feedback delays the 
growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely. 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) are among the most vio¬ 
lent explosions in the Universe. They are classified into 
two populations by the duration of prompt emission Tgo: 
short GRBs (SGRBs) wit h Tgn < 2s and long G RBs 
(LGRBs) with Tgo > 2 s (jKouvelioton et all 1199, Sll . A 
widely accepted model of LGRBs is the collapsar sce¬ 
nario, in which a black hole (BH) and accretion disk sys¬ 
tem is formed after the stellar collapse, and launches a 
relativistic jet th at breaks out the progenitor star, pro- 
ducing a GRB (IWooslevI 119931 IMacEadven fc WooslevI 
[TMl). Theoretical models identify the typical dura¬ 
tion of a LGRB (Tgo ^ 30 s) as the free-fall time of 
the envelope, or the sound-crossing time of the shocked 
envel ope, of a Wolf-Rayet (WR) star (|Mizuta fc lokal 
1201311 . Erom observation we know that at least some 
LGRB s are accompanied by b r oad-lined Ic supernova e 
(SNe; IWooslev fc BloomI 120061 iHiorth fc BloomI 1201211 . 
This suggests a tight connection between GRBs and pro¬ 
genitors with stripped envelope like WR stars. 

Recently, some LGRBs have been discovered to show 
ultra-long duration of the prompt emission with St.y ^ 
10^ s. These have b een named ultra-long gamma- 
ray b ursts (ULGRBs) (|Gendre et al.l 120131 ILevan et al.l 
1201411 . The ultra-long duration was first pr edicted in the 
context of Population HI ( PopHI) GRBs (iSuwa fc lokal 
120111 iNagakura et al.l[2012ll , and subsequent studies sug¬ 
gested that a metal-poor blue supergiant (BSG) collap- 
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sar rather than t he WR one is more favorable to ex- 
plain such bursts (IGendre et al.ll2013l iKashivama et al.l 
120131 iNakauchi et al.l 1201311 . Since metal-poor stars 
may suffer from little mass loss, metal-poor BSG stars 
would keep mas sive hydrogen envelo pes until the pre¬ 
collapse phase (|Wooslev et al.l 1200211 . Therefore, the 
accretion of the massive hydrogen envel ope can lead 
to th e long-lasting cen t ral en g ine activity (ISuwa fc lokd 
20111 INagakura et al.l 120121 iQuataert fc Kasenl 120121 
Wooslev fc Hegeii l2012tlNakauchi et al.l 1201311 . Since 
metal-poor stars are considered to be the dominant pop¬ 
ulation in the high-z Universe, ULGRBs might be a dom- 
inant population of GRBs in the high-z Universe (e.g., 
Ide Souza et al.ll^OllH . 

On the other hand, the existence of BHs of mass 
^ 10® M|^ at z > 6 (|Eanl 120061 iMortlock et al.l 120111 
iWu et ani2015ll is a great mystery in astrophysics be¬ 
cause the accretion time is not enough to grow the 
BHs from the stellar mass BHs of mass ^ 10 Mq (e.g., 
iHaimanl 1201311 . Many attempts are made to solve this 
problem. One possible solution may be the forma¬ 
tion from PopHI stars. Recent numerical simulations 
suggest th at the mass of PopHI stars reaches up to 
10®“® Mq (jHirano et al.l 1201^ iSusa et al.l 1201411 . BHs 
born from these PopHI stars can barely grow up to su¬ 
permassive BHs of mass 10® Mq at z > 6, if succes¬ 
sive high mass accretions are maintained. Feedback ef¬ 
fects from accreting BHs, how ever, decrease the accre¬ 
tion rate (lAlvarez et al.l [200^ . Thus it seems difficult 
for seed BHs from PopHI stars to form the observed su¬ 
permassive BHs without resorting to super-Eddington 
accretion via such as Bondi and cold acc retion (e.g., 
iVolonteri fc Reej 120051 iGaspari et al.l l2013l l or efficient 
mergers (e.g., iMadau fc Reesll2001ll . 

An attractive alternative might be the formation of 
supermassive stars (SMSs) of mass ~ 10® M©. SMSs 
may be formed in the high temperature region irradiated 
by the strong ultra-violet radiation from nearby galaxies 
(I0mukaill20011 IB romm fc LoebI 120031 iShang et al.l 120101 

Latif et al.l 120131 llnavoshi et al.l 12014 [Johnson et al.l 
201411 . When SMSs end their life due to the ex¬ 
haustion of their nu clear fuel or the general relativis- 
tic (GR) instability (|Chandrasekhad[l964 lOseiO 119661 
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iShapiro fc Teukolskvl 119831 : iShibata &: Shapir^ I2002D . 
they collapse to massive BHs of mass ^ 10® Mq, so-called 
direct collapse BHs (DCBHs). If these DCBHs are seeds, 
the accretion time might be enough to grow the BHs to 
- 10 ^ M 0 at z > 6 . 

Various violent phenomena from the gravitational col¬ 
lapse of SMS s are expected such as the energetic neu¬ 
trino bursts (jFrver fc Heererl 1201 111 and SNe with very 
huge explosion energy of ^ 10®® erg ll Johnson et al.l[2013l : 

I Whalen et al.l I2013t IChen et al.l 1201411 . In this paper, 
we consider ULGRBs from SMSs as another possibil¬ 
ity. If SMSs evolve without mass ejection and collapse 
to form BH-disk systems, like collapsars, we can ex¬ 
pect the launch of the relativistic jet similar to LGRBt@. 
Since SMSs are larger in radius than PopHI first stars 
of mass 10-1000 Mq or BSGs, we expect that the dura¬ 
tion of ULGRBs is even longer than that of the observed 
ULGRBs. The detection of such ULGRBs would en¬ 
able us to observe the very moment of the birth of first 
quasars to probe the high -2 Universe. Such energetic ex¬ 
plosions might prevent the subsequent gas accretion and 
the growth of the DCBHs. Thus, it is worth evaluating 
the details of relativistic jet explosions from supermas- 
sive collapsars and their observational signatures. 

This paper is organized as follows: In we show the 
pre-collapse stellar models of SMSs, which evolve from 
zero age main sequence (ZAMS) stars. In iJSl we de¬ 
scribe the method to calculate the jet dynamics in the 
SMS envelope. In SJU we show that the jet can break out 
the SMS. We also discuss the observational signatures 
and the detectability of GRBs from supermassive collap¬ 
sars. In lj5l we discuss the other progenitor model in 
which the SMS is accreting mass and collapses through 
the GR instability. Then we estimate the event rate of 
GRBs from SMSs, and discuss the effects of GRBs on 
their environment. A summary and our conclusions are 
given in SJ5] Throughout this paper, we consider the 
ACDM cosmology and adopt the cosmological param¬ 
eters as : i Jn = 67.8 km s~® Mpc~®, Dm = 0.308 and 
Da = 0.692 (jPlanck Gollaboration et ^120151) . 


2 . PROGENITOR MODEL 


Theoretical studies of the primordial star forma¬ 
tion have suggested that SMSs with > 10® Mq 

can be formed under the hot environments of 


first galaxies (lOmukail 120011: IBronim & LoebI 

2003 

Shansf et al.ll2010l: Latif et al.ll2013l: llnavoshi et al.l 

2014 

Johnson et al.l l2014fl. Whereas molecular hvdroEen is 


the primary coolant in primordial star-forming clouds, 
its formation is prevented under strong ultra-violet ra¬ 
diation from nearby galaxies, so that the temperature is 
kept at 10^ K in such a cloud due to the hydrogen 
atomic cooling. The accretion rate onto a protostar M 
can be evaluated by dividing the Jeans mass Mj of the 


® When a SMS collapse to a BH, there is some possibility that 
a qua si-star, which is power ed by the central BH accretion, could 
form UBegelman et alJ]2008ll . In this work, however, we consider 
that t he central BH launches a jet as o rdinary GRBs. Previous 
works IIBarkovI 120101: fCzernv et al.ll2012l) did not consider the jet 
propagation in the SMS, so that the observed quantities such as 
the luminosity and duration are unable to obtain. 


star-forming cloud by the free-fall time iff: 

/ T \ 

M ^ Mj/fff ^ cl/G ^ 0.51 Mo yr-i, (1) 

where Cs and T are the sound velocity and the temper¬ 
ature of the cloud, and G is the gravitational constant, 
respectively. We also assume that the cloud is composed 
of only hydrogen, for simplicity. Thus, the mass accretion 
rate onto the protostar formed in the central region of the 
star-forming cloud can be as high as ^ 0 . 1-1 Mq yr“^ un¬ 
der such hot environments. The protostar can grow up 
to a SMS with > 10® M 0 within its lifetime of ~ 1 Myr 
through such a high mass accretion rate. 

When the mass accretion onto the protostar stops be¬ 
cause of radiation feedback, or for other reasons, the pro¬ 
tostar contracts onto the ZAMS in a Kelvin-Helmholtz 
timescale. Since the SMSs are almost fully convective 
and the Kelvin-Helmholtz timescale is short compared to 
the lifetime of the star, the resulting final stellar struc¬ 
ture should be the same as long as accretion stops not 
too close to hydrogen depletion. After the SMS runs 
through its nuclear bur ning phase^lB, it collapses to a BH 
similar to a massive staiFTWhen the SMS collapses to a 
BH, it may produce a GRB by launching a relativistic jet 
provided there is enough angular momentum. Thus, we 
adopt this pre-collapse SMS for our progenitor model. 

In this paper, we focus on a supermassive progeni¬ 
tor with a ZAMS mass of 10® M 0 , but for comparison 
also con sider a progenitor wh ich has the ZAMS mass of 
10^ M 0 (|Erver fc Ilegeil[2011ll . We call the former model 
as “1E5 model” and the latter model as “IE4 model”, 
respectively. The density profiles of these pre-collapse 
SMSs are shown in Eigure[ll We can see that they have 
very large radii of i?* ^ 10 ^^ cm, which are as large 
as those of present-day RSGs (green curve in Eig. [1]). 
Note that the envelopes of the SMS models have steeper 
density profiles {p oc r“®) than that of the RSG model 
(p oc This is because radiation pressure domi¬ 

nates in the SMS envelopes while gas pressure dominates 
in the RSG envelope (see Section lOl) . 

It is not trivial whether the relativistic jet can break 
out the envelope of the supermassive progenitors suc- 
cessful ly since it can not break out the large envelope of 
RSG s ( Matzneill2003l:lSuwa fc IoO201IllNakauchi et al.l 
I 2 OI 2 II . ISuwa fc lokal (i2011h : lNagakura et al.l ( 201211 stud- 
ied relativistic jet explosions from massive PopHI stars 
(M, ~ 10 ® M 0 ), and found that the jet breakout is pos¬ 
sible despite the large radius of i?* ^ 10®® cm. They at¬ 
tributed this to the long duration of the mass accretion of 
the massive envelope. In this paper, we consider super¬ 
massive PopHI stars which have M* ~ 10"®-10® M 0 and 
Rt ^ 10 ®^ cm to determine whether the jet can break 
out of the envelope. These stars have 10 times larger 
ra dii compared t o the massive PopHI s t ars co nsidered 
by ISuwa fc lokal (|2011ll : iNaeakura et al.l (|2012ll so that 
the jet might not break out the envelope. 

3. NUMERIGAL METHOD 
3.1. Jet Model 

® The advanced phases may be accelerated as part of the collapse. 

Even a p hase with a hot pr oto-BH due to neutrino trapping 
as observed bv IPrver et al.l 11200 11 1 may nor occur. 
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Figure 1. Density profiles of SMSs in their pre-collapse phase. 
The red curve and blue curve correspond to the 1E5 model and 
the 1E4 model, respectively. These profiles are obtained by cal¬ 
culating the metal-free 10^ Mq and 10^ Mq ZAMS stars until the 
nuclear fuels are exhausted in the center. The green curve shows 
the d ensity profile of a RS G (pre-supernova structure of Model sl5 
from [WooslW et al.l 12003^ . These profiles show that SMSs have 
similar radii with a RSG but steeper density profile than that of a 
RSG. The gray dash-dotted lines show the slope of the power-law 
density profiles of p oc r~^ and p oc as references. 


We adopt a simple jet model in which the jet luminosity 
varies with ti me and depends on t he mass accretion rate 
onto the BH (|Suwa fc ^okall2011^ . The SMS begins to 
collapse, when its nuclear fuel is exhausted in the center, 
and forms a massive BH. We evaluate the mass accretion 
rate onto the central BH using the free-fall timescale as 
typical timescale. For each mass shell at radius r and 
mass coordinate Mr, the free-fall timescale can be calcu¬ 
lated from 



where p = Mr/^r^ is the mean density within r. The 
mass accre tion rate is calculated by (iWooslev fc HeeeJ 


M = 


dMr dMr/dr 2Mr 
dtff 


dts/dr ts{r) \p-p 


( 3 ) 


A relativistic jet will be launched from the cen¬ 
tral engine, which is composed of a BH and an 
accretion disk. The mechanisms of launching rel¬ 
ativistic outflows are still uncertain, though several 
processes have been proposed such as MHD mech¬ 
anism and neutrino-antineutrino anni h ilation mecha¬ 
nism ( e.g., IBlandford fc ZnaiekI Il977t iPopham et al.l 
1199911 . ISuwa fc lokal ( ~llll studied the jet-driven ex¬ 
plosions from massive PopHI stars and found that the 
MHD mechanism works long enough for the jet to break 
out the envelope, whereas the neutrino mechanism does 
not. Therefore, we assume the MHD mechanism here¬ 
after. In the MHD p rocess, the jet luminosity ca n be 
modeled according to ([Komissarov fc Barkovll20lOl 


Lj(t) = ?7jM(t)c2, (4) 

where rjj is the efficiency parameter. For simplicity we 
assume that rjj is constant. We can take into account the 


effect of the progenitor rotation by adjusting the parame¬ 
ter value ? 7 j. When a sufficiently fast-rotating progenitor 
collapses, it can form a BH and accretion disk system. 
The remainder of the star then falls onto the central 
BH in the accretion time face which can be related to 
the free-fall ti me as face oc wh ere a is the viscos- 
ity p arameter (iShakura fc Sunvae^ 119731 : iKumar et al.l 
120081 1. Thus, the factor a can be absorbed into the ef¬ 
fective value of 77j. We use pj = 6.2 x 10“^, which is 
calibrated to reproduce the total jet energy of typical 
LGRBs of Etnt = f mMc^dt = 1 0^^ erg for a WR pro¬ 
genitor model (ISuwa fc Iokall201lll . 

Throughout this section, the time t is measured in the 
central engine frame. We set t = 0 as the time when the 
black hole is formed. Since the formation mechanism of 
a relativistic jet is still under debates, the time of jet for¬ 
mation, tin, is also uncertain and should be a parameter 
of this study. Following the previous studies, we assume 
that the jet is formed when the mass of BH reaches 3 Mq: 
tin = iff(Dn), where the enclos ed mass within the radius 
Tin is 3M0, i.e., = 3M0 (jSuwa fc Iokall20lTl l. The 

choice of tin (or however, has little influence on 

the jet dynamics in the envelope as long as tin is much 
smaller than the jet break out time. In fact, we find that 
the jet breakout time is within a factor of two, when the 
jet formation time is changed as = 30 M 0 , 300 M 0 , 
and 3,000 M 0 in the 1E5 model. 


3.2. Jet Propagation in the SMS Envelope 

Here, we describe the jet dynamics in the progen¬ 
itor e n velope followi n g the prescription in Matznen 
(1200311: ISuwa fc loka (1201 If) : iBromberg et al. ( 201111 : 
iNakauchi et al.l ( 2012 1. A relativistic jet launched from 
the central engine collides with the stellar matter and 
forms the shocked region at the jet head. The jet head is 
composed of the forward shock which sweeps the stellar 
matter and the reverse shock which decelerates the jet 
matter. Both shocked matter are divided by the con¬ 
tact discontinuity. Then the velocity of the jet head can 
be calculated from the pressure balance at the contact 
discontinuity as 

- Pirf + Pj = PnC^hJdlPl + Pa, (5) 


where, p, h, and P represent the mass density, specific 
enthalpy, and pressure measured in the fluid rest frame, 
respectively, and j5 and F = (1 — are the veloc¬ 

ity normalized by the speed of light c and the Lorentz 
factor, respectively. The subscripts “h”, “j”, and “a” in¬ 
dicate that the quantity is measured in the rest frame of 
the jet head, the jet, and the ambient stellar medium, 
respectively. Since the stellar medium is composed of 
the non-relativistic matter, we can neglect its pressure 
against the rest mass energy density, i.e., Pa ^ PaC^, 
and this also leads to ha — 1. In the left hand side of 
Equation ([5|), we can also neglect the jet pressure, since 
the jet is ultra-relativistic (/3j ~ 1 and Fj 1). 

Using the above approximations, the velocity of the jet 
head is given by 


/3h 


1 


l + L- 


( 6 ) 


where L = p-^h-^T^/p^, is a parameter which determines 
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the jet dynamics. Using the collimation-corrected jet lu¬ 
minosity given by 




(7) 


L is given by 


L ~ 


PaC^Sh’ 


( 8 ) 


where Eh is the cross section of the jet head. Thus, L 
is the ratio of the luminosity of the jet to that of the 
ambient stellar medium rest energy flux. Here, we as¬ 
sume that the opening angle of the jet is constant with 
9 = 5°. This is the typical val ue obtained from the af¬ 
terglow observations of LGRBs (iFrail et al.ll2001h . Then 
the cross section can be given by Eh = 7r(rh0)^, where 
rh{t) = f cPhdt is the radius of the jet head. 

As long as the velocity of the jet head is non- 
relativistic (/3h < 1), the shocked jet head can expand 
sideways to form a cocoon structure surrounding the 
jet. Since the temperature of the cocoon is high, it 
is radiation-pressure-dominated. As long as the sound 
crossing time in the cocoon is shorter than the dynamical 
time of the jet head, we can neglect the inner structure 
of the cocoon and assume that the cocoon is uniform. 
Hereafter, we consider the one-zone model for the lateral 
expansion of the cocoon. The cocoon is overpressured 
with respect to the ambient stellar medium so that it ex¬ 
pands laterally. By considering pressure balance at the 
surface of the cocoon, the l ateral expansion velocity of 
the cocoon is calculated by (jBegelman fc Cio^ll989l) 


/3c = 



(9) 



Figure 2. Time evolution of the velocities of the jet head and 
the lateral expansion of the cocoon surface. Red, blue and green 
lines correspond to the 1E5 model, the 1E4 model and the RSG, 
respectively. The horizontal axis shows the time from jet formation. 
The vertical axis shows the velocity divided by the speed of light. 
The solid and dashed curves correspond to the velocities of the 
jet head and the cocoon surface /3c, respectively. For the 1E5 
model and the 1E4 model, the jet head velocity is always larger 
than that of the cocoon edge so that the supermassive collapsar 
jet can break out the envelope. For t > 500 s and t > 2,000 s, 
the jet head is accelerated drastically in the 1E5 and 1E4 models, 
respectively, because the envelope density decreases more steeply 
than oc r“^. On the other hand, in the RSG, the jet head velocity 
is overtaken by that of the cocoon at ~ 3,000 s. 


Substituting Equations (HJ and m and the definition 
of Eh into Equation ([6]), the jet head velocity is given by 


/ 7rc6)Va(r-h)rg \ ^ 

V VjMit) ) 


( 12 ) 


where Pc is the pressure in the cocoon and Pa(?'h) = 
Mj.^/(47rrj(/3) is the mean density of the progenitor star 
within the radius rh. 

Since the cocoon matter is radiation-pressure- 
dominated, the pressnre is given by Pc = Ec/iVc, where 
Ec and 14 represent the total energy and volume in the 
cocoon, respectively. Becanse the cocoon energy is snp- 
plied from the jet head, it is given by 

Ecit)=T]cJ^ ( 10 ) 

where rjc indicates the fraction of the jet luminosity 
streaming into the cocoon. Throughout the paper, we set 
rjc = 1, since the velocity of the jet head is non-relativistic 
for most of the time within the progenitor envelope. In 
Equation (|TTHl . the upper limit of the integral indicates 
that at t, the jet head receives the luminosity which is 
produced at f — rh{t)/c at the central engine. For the 
cocoon volume, we assume that the cocoon has a coni¬ 
cal shape with the height of rh(t) and the base radius of 
rc(<), so that it is given by 

Veit) = Krl{t)rh{t), ( 11 ) 

where rdt) is the lateral distance of the cocoon surface 
from the jet axis given by rc(t) = J cfdcdt. 


Substitnting Equations cni) and (ED into (|9]), the lateral 
expansion velocity of the cocoon is given by 


4. ULGRBS FROM SUPERMASSIVE COLLAPSARS 

4.1. Successful Breakout of Supermassive Collapsar 

Jets 

In Figured! we show the time evolution of the velocities 
of the jet head and the lateral expansion of the cocoon 
surface with solid and dashed curves, respectively. Each 
color corresponds to the 1E5 model (red), the 1E4 model 
(blue), and the RSG model (green), respectively. They 
are calculated from Equations (fT^ and (fT^ . The hor¬ 
izontal axis shows the time from jet formation: t — <in. 
The vertical axis shows the velocity divided by the speed 
of light. 

We can see that for the models 1E4 and 1E5, the ve¬ 
locity of the jet head is always larger than that of the 
cocoon surface. On the other hand, for the RSG model, 
the velocity of the cocoon exceeds that of the jet head at 
> 3,000 s. In the latter case, the radius of the jet head 
is comparable to the lateral size of the cocoon surface so 
that they can reach the stellar surface almost at the same 
time. This looks like a spherical explosion rather than a 


° Li{t')dt' 

3{Mry^ - 


J./Z 


(13) 
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Figure 3. Time evolution of the mass accretion rate onto the BH. 
The horizontal axis is the time from jet formation. The red and 
the blue lines corresponds to the accretion rate obtained from the 
1E5 model and the 1E4 model, respectively. Whole in each color, 
dashed and solid lines correspond to the time before and after the 
jet breakout, respectively. The shaded region in each color means 
that the jet head propagates in the progenitor’s interior and the jet 
is “hidden” from the observers. The accretion rate obtained from 
the WR model corresponding to LGRBs is shown by the gray curve 
as a reference. We can see that supermassive collapsar jets can lead 
to ULGRBs with the duration of Tgo ~ 10"* s up to > 10® s. Note 
that the light curves decline steeply after the breakout, so that the 
observed duration Tgo is approximately equal to the breakout time. 


collimated explosion. On the other hand, in the former 
case, the collimated jet breaks out the progenitor surface 
since the jet head reaches the surface much earlier than 
the cocoon. 

In Figure|3l we show the time evolution of the mass ac¬ 
cretion rate onto the BH. The horizontal axis represents 
the time from jet formation. The red and blue curves 
correspond to the accretion rate obtained from the 1E5 
model and the 1E4 model, respectively. While on each 
curve, dashed and solid regions correspond to the time 
before and after the jet breakout, respectively. We also 
show the accretion rate obtained from the WR model 
(i.e., LGRB) with the gray curve as a reference. We can 
see that supermassive collapsar jets can lead to ULGRBs 
with the duration of > 10^ s, owing to the accretion of 
the massive envelope. 

While the SMS models have as large radii as 
present-day RSGs, we find that the relativistic jet 
can break out the progenitor envelope successfully. 
This can be attributed to the difference in the slope 
of the density profile. At the outer envelope of 
a polytropic star with the polytropic index n, the 
density profile can be approximated as P air) oc 
{R*/r — 1)" ~ r~^ (|Matziier fc McK^ 119991 1. Since 
SMSs have radiation-pressure-dominated convective en- 
vel opes, the index can be a pproximated as n = 
3 (iShapiro fc Teukolskvl 119831 : iKippenhahn fc WeigertI 
Il990li . On the other hand, RSGs have gas-pressure- 
dominated convective envelopes so that the index is 
n = 1.5. Thus, the SMS models have steeper density 
profiles than the RSG models. 

When the velocity of the jet head is non-relativistic, L 
is less than unity, so that the velocity can be approxi¬ 


mately evaluated as 

/3h(t) oc (14) 

Furthermore, by using that M is approximated as 

M(t) oc t~3~ in t he density profile of Pa{r) oc 
(|Snwa fc lokal 120111 1 ■ and that rh is evaluated roughly 
as Th oc /3ht, we obtain 

/3h(t) oc rh(t)5^, (15) 

from Equation ([14]). Therefore the jet head is acceler¬ 
ated in the outer region of the SMS envelopes where the 
density profile decreases more steeply than oc r~^. On 
the other hand, the jet head is decelerated in the RSG 
envelopes. From Figured we can see that the jet heads 
are accelerated drastically at f > 500 s and t > 2,000 s 
for the 1E5 and 1E4 models, respectively, when the jet 
heads enter the regions where the envelope density de¬ 
creases more steeply than oc r~^. On the other hand, 
for RSGs, the velocity of the jet head is decelerated in 
the envelope so that it takes much longer time to break 
out the envelope. In this case, the lateral size of the co¬ 
coon becomes comparable to the radius of the jet head. 
This is a spherical explosion but not a collimated GRB. 
Th us, RSGs canno t be the progen i tor of LGRBs as shown 
bv iMatzned (|2003l l: ISuwa fc lo^ (j20lH ): lNakanchi et al.l 

(l20m 

It is possible that a disk wind flows out isotropically 
from the accretion disk and changes the envelope struc¬ 
ture. In this case, the disk wind and the deformed en¬ 
velope may affect the jet propagation. Unless the wind 
velocity is larger than the jet head and cocoon veloc¬ 
ity, however, this effect can be ignored. The wind veloc¬ 
ity is evaluated from the wind energy as ^ 
where is the mass of the wind component. The wind 
energy is also given by ~ where is the 

efficiency parameter, as defined for the jet luminosity in 
Equation ([I]). On the other hand, the cocoon velocity 
is evaluated by E^ ^ Mcv1/2 J rjjMc^dt. From these 
equations, we obtain the ratio of the cocoon velocity to 
the wind velocity as Vc/v^ ^ . In our 

work, the jet efficiency parameter is ? 7 j ^ 6 x 10“^. As 
the wind flows out isotropically while the cocoon expands 
around the jet, their masses are related as Me ^ /2, 

where 0 0.1 is the jet opening angle. Then, we have 

Vc/vw ~ 0.377w ' ■ Thus, when the disk wind flows out 
very efficiently [rj^ > 0.1), we should consider the effect 
of the isotropic wind outflow on the envelope structure 
and the jet propagation. The recent numerical simu¬ 
lations of super-Eddi ngton accretion d isks suggest the 
value of 77w ~ 0.01 (iJiang et al.l l2014lf'^ . which yields 
Vc 3i;w It should be noted that the accretion rates in 
their simulations are M ~ lOOLEdd, while in our situa¬ 
tion, the mass accretion rates amount to M > lO^^LEdd- 
We need numerical calculations in order to study the disk 
wind in our case. 

4.2. Prompt emission 

ISadowski et al.l II2014I ) also gives the efficiency ~ 0.3, although 
this value includes the efficiency from the accretion energy to the 
radiation and magnetic energy. 
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Once the relativistic jet breaks out the progenitor’s en¬ 
velope, it can contribute to the prompt high-energy emis¬ 
sion, like LGRBs. Here, we evaluate the observational 
signatures and the detectability of the prompt emission 
from the supermassive collapsar jets. Since the emission 
mechanisms of the prompt emissi on of GRBs are s till un - 
der debates, however, following iNakauchi et al.l (I2012I1 . 
we evaluate them by applying simple empirical relations 
to supermassive collapsar jets. 

First of all, we describe our model for the prompt emis¬ 
sion. We assume that once the jet breaks out the surface 
at tb) the relativistic jet can contribute to the gamma-ray 
emission by using a fraction of its energy. Thus, the 
collimation-corrected gamma-ray luminosity is given by 
L^{t) = e^Lj{t). We also assume that the high energy 
emission lasts until all the matter in the envelope has 
accreted onto the central BH at Hence, 

the duration of the prompt emission can be evaluated 
by tg,* — th- Hereafter, we adopt 9 = 5°, = 0.1 

and T]j = 6.2 x 10“"^ as our fiducial values. Then the 
isotropic luminosity of the prompt emission is given by 
L-y,iso{t) = . In this model, the isotropic 

radiated energy of the prompt emission E^^iso and the 
peak luminosity Lp are estimated by 


rifi,- 


£^7,iso— ^2 / L^{t)dt — 


'th 


6>2 


rUi, 


'th 




= 2.9 X 10'’ 


52/ Jtb 




IMm 


ergs. 


(16) 


and 


Z 9 • 

-Zvp = Z/^pso(^ — ^b) — ^ “ ^b) 


= 2.9 X 10^2 

V IMq s-1 


ergs s 


-1 


(17) 


As long as the luminosity is proportional to the mass 
accretion rate, the luminosity decreases monotonically 
after the breakout. Therefore, the luminosity peaks at 
the breakout. 

Next, we evaluate the spectral peak energy of the 
promp t emission in the cen tral-engine frame Ep. Fol¬ 
lowing INakauchi et al.l (l2012tl , we assume that either one 
of the two empirical correlations of LGRBs holds in su¬ 
permassive collapsar jets: the ifn-Tn correlation or the 
E^-E ^ isn correlation (jYonetoku et al.l[2004l: lAmati et al.l 
l2002fl . If the Ep-Lp correlation holds for the burst, then 
the spectral peak energy can be evaluated from the peak 
luminosity Lp, using the correlation 


L„ 


1052 gj.g g 


~ 2 X 10 


-5 




1 keV 


2.0 


(18) 


as Ep = 5.6 X 10^ and 1.6 x lO^keV for the 1E5 model 
and the 1E4 model, respectively. On the other hand, if 
the Ep-Ejpso correlation holds for the burst, Ep can be 
evaluated from the isotropic radiated energy E^^so, using 
the correlation 


Ep 

1 keV 



E • 

-^7,ISO 

1052 


0.57 


(19) 


as Ep = 2.6 X 10'^ and 1.4 x 10"^keV for the 1E5 model 


Table 1 

Observational Characteristics of the Prompt Emission at 2 = 15 


Progenitor Model 

IL 5 

1L4 

Ejpso [erg] 

Lp [ergs-l] 

2.5 X 10”'’ 

6.2 X 10^2 

8.4 X lO”” 

5.1 X 10®* 

[keV] 

3.5 X 10 

1.0 X 10 

Eg*’® [keV] 

1.6 X lO’* 

8.6 X 10^ 


Notes. in line 4 and in line 5 show the peak energy 

of the spectrum predicted by the empirical relations lO and mil, 
respectively. 

and the 1E4 model, respectively. We summarize the ob¬ 
servational signatures of the prompt emission from the 
supermassive collapsar jets in Table [TJ where the redshift 
of the burst is set as 2 ; = 15. In lines 4 and 5, the spec¬ 
tral peak energy in the observer frame Lp'’’® is given by 
^obs _ Ep/{1 + z). From Table [U we find that the to¬ 
tal energy is much larger than that of LGRBs, while the 
peak luminosity is comparable to them. The accretion 
time of SMS is much longer than that of LGRB so that 
SMS releases much larger amount of energy than WR 
collapsars although the luminosity is similar. 

Finally, we discuss the detectability of the prompt 
emission from the supermassive collapsar jet with detec¬ 
tors like the B urst Alert Telescope (B AT) on board the 
Swift satellite (|Barthelmv et akl 120051) . BAT covers the 
energy range from Lmin = 15keV to Lmax = 150 keV. 
The energy flux detected by BAT is given by 

EN{E)dE 

/sig(7,obs) = ^b°i(b) ffEN{E)dE ’ 

where tj = t — tg, 7,obs = (1 + z)t^, Ef{E) and Lboi(7) 
are the time from the breakout, the time in the ob¬ 
server frame, the photon number spectrum and the bolo- 
metric flux, respectively. Empirically, we assume that 
N(E ) is represented by the Band function (iBand et al.l 
1199, *!() with the spectral indices of a = —1 and /3 = 
—2.3 (IKaneko et al.l[200^ . The bolometric flux is given 
by 

Lboi(t 7 ,obs) = erg s“^ cm“^, (21) 

where di^[z) is the luminosity distance. 

In Figure 01 we show the light curves of the prompt 
emission in the case of the Ep-Lp correlation. The red 
and blue curves correspond to the 1E5 and the 1E4 mod¬ 
els, respectively. We set the redshifts of the bursts as 
z = 10 (solid), 15 (dash-dotted), 20 (dotted), respec¬ 
tively. The gray dotted lines show the BAT sensitivities 
/sen(Atobs) with the integration times of Atobs = 1 s, 
10^ s, and 10^ s, from top to bottom. If the burst enters 
the field of view of the BAT at some time t^^obs; and the 
signal flux is larger than /sen(Atobs), then it can be ob¬ 
served by BAT up to 7,obs + Atobs- We can see that the 
burst can be detectable up to z = 20 for Atobs = 10^ s 
in the 1E5 model. We can also see that the burst can 
be detectable up to z = 20 for Atobs = 10^ s in the 1E4 
model. 

In Figure 01 we show the light curves of the burst in 
the case of the Ep-Ej^so correlation. We can see that 
the observed flux is smaller by an order of magnitude 
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Figure 4. Light curves of the prompt emission of the ULGRBs 
from supermassive collapsars. The flux is calculated at the Swift 
BAT energy range (15-150 keV), assuming that the Ep-Lp correla¬ 
tion holds. The red and blue curves correspond to the 1E5 and the 
1E4 models, respectively. The redshifts of the bursts are z = 10 
(solid), 15 (dash-dotted), 20 (dotted), respectively. The gray dotted 
lines show the BAT sensitivities /sen(Atoi,s) with the integration 
times of At^ns = 1; 10^ i 10^ from up to bottom. 



Figure 5. The same as Figure |4] but the flux is calculated as¬ 
suming that the Ep-E-ypao correlation holds. 


than that of the above case. This is because the Ep- 
£' 7 ,iso correlation leads to ~ 1 MeV, which is out of 
the BAT energy range. Nonetheless, we can see that the 
burst can be detectable up to z = 20 for Atobs = lO'^s 
in the 1E5 model. On the other hand, in the 1E4 model, 
longer integration times (Atobs > lO'^s) are needed to 
detect the burst at such a high redshift of 20 . 

Thus, we conclude that the prompt emission from su- 
perinassive collapsar jets can be detectable as ULGRBs 
up to z ^ 20 even with the current detectors like BAT. 

5. DISCUSSION 
5.1. Accreting supermassive stars 

If the high mass accretion rate of 0.1-1 Mq yr“^ con¬ 
tinues, the SMS mass will eventually reach a few times 
^ 10^ Mq and the GR instability sets in to lead the for¬ 


mation of a BH and accretion disk system H- Such a 
high mass accretion rate onto protostars is obtained if 
the temperature is high as shown in E quation ilL and 
second paragraph of Section[2j Recently. iHosokawa et al.l 
(1201,11) calculated the evolution of a SMS under such a 
high mass accretion rate ^ 0 . 1-1 M 0 yr“^ from the pro- 
tostar phase to the final mass of M* = 10® Mq. They 
found that the SMS evolves with a very large envelope of 
i?* ~ 10®® cm. Under the rapid mass accretion, the stel¬ 
lar envelope has a large opacity dominated by H~ ions, 
and absorbs heat released by the contraction of the stel¬ 
lar inner region. Then, the st ellar envelope expand s as 
t he protostar accretes in atter (jHosokawa et al.l[2012D . 

IHosokawa et al.l (j 20 l, 1 f l. however, were not able to com¬ 
pute the SMS evolution beyond 10® M 0 , because their 
stellar evolution code is suffered from the numerical dif¬ 
ficulties. The reason is not clear. If the accretion is 
stopped before the GR instability sets in, the final pre¬ 
collapse model is similar to the 1E5 model so that the 
ULGRB of the SMS is expected as discussed in the pre¬ 
vious section. If the accretion continues, the SMS will 
enter the GR instability region. In such a case, we can 
ass ume that the SM S has s imilar density profile obtained 
by IHosokawa et al.l (j20l,1[) when it begins to collapse 
through the GR instability. The structure of the SMS 
envelope would not change so much as long as the high 
mass accretion rate is kept until the SMS obtains the crit¬ 
ical mass. Therefore, our assumption may be justified for 
the envelope, which is important for the propagation of 
jet heads as we saw in the previous section. 

In Eigureini we show the density profile of the accret¬ 
ing SMS when its mass reac hes M^, = 10® Mf^ by using 
the magenta curve. While in IHosokawa et all (I2013D . the 
density profile is given as a function of the mass coor¬ 
dinate Mr (in Eig. 3 of their paper), we show it as a 
function of radius r by integrating the mass conservation 
equation dr/dAfr = l/ 47 rr^p(A'/r). As shown in Eigure 6 , 
the accreting SMS is about 10 times larger than a RSG. 
Hence one may expect that the jet can not break out the 
surface of the star as RSGs. However, we find this pro¬ 
genitor also has a steep density profile in the radiation- 
pressure-dominated env elop e. Therefore, according to 
the argument in Section [TTJ we expect that a relativistic 
jet is accelerated to break out the surface of the progen¬ 
itor star. We calculate the jet propagation in the same 
way described in Section [3l and find that the jet actu¬ 
ally breaks out the envelope. Thus, we expect that even 
if the SMS collapses in the accreting protostar phase, it 
can produce an energetic explosion. We also calculate 
light curves of the prompt emission, assuming that the 
Ep-]^ and the Ep-E-y^iso correlations hold. In Figures [7] 
andIHl we show the light curves. ULGRBs from accreting 
SMSs are dimmer than the GRBs from the 1E5 model. 
This is because the accreting SMSs are larger than SMSs 
of the 1E5 model. It takes more time for a jet to break 
out the larger envelope. Thus, when the jet breaks out 

If a SMS rotates , it is stabilized against th e GR 
instability (e.g., IFow leij 19661 : IBisnovatvi-Kogan et all 119671 : 
IBaumgarte fc ShaDirolll99E ). However, even if the SMS rotates, 
it will not acquire more mass than ~ 10® Mq for the highest mass 
accretion rate of 1 Mqs“^. Eventually, it exhausts the nuclear 
fuel in ~ 1 Myr and collapses to form the central BH. Actually, 
the calculation of lHosokawa et aH ll201.'il shows that the hydrogen 
burning occurs in the core. 


























































Figure 6. Density profile of an accreting SMS with mass of 
10 ^Mq {magenta curve). This progenitor also has the steep den¬ 
sity profile of the envelope (with slope of ~ 3), where a relativistic 
jet is not decelerated. Note that this progenitor is not yet in the 
pre-collpase phase, but in the accreting protostar phase. We can 
expect that when the GR instability sets in, however, the pre¬ 
collapse progenitor has the similar density profile especially in the 
envelope as long as the high mass accretion rate is kept. With the 
red and green curves^ we also show the density profiles of the 1E5 
model and RSG, respectively, for references as in Figure^ 



Figure 7. The same as Figure [4] but the flux is calculated from 
the accreting supermassive collapsar {magenta curves) assuming 
that the E-p-L-p correlation holds. We also show the light curves of 
the 1E5 model for references. 

the envelope, the accretion rate onto the central BH is 
decreased. The observation time to detect ULGRBs from 
accreting SMSs is longer than that of the non-accreting 
SMSs in Figures m and [S] 

5.2. Event Rate 

We briefly discuss the detection rate of the ULGRBs 
from supermassive collapsars. For a given observation 
time Atobs, the cumulative number of ULGRBs AN{z) 
which have redshifts less than z can be calculated from 

AN{z)= [ '^GRB{z')^7Tcr{z')‘^ 

Jo 

where 'I'grb is the intrinsic event rate of ULGRBs and 
r(z) is the comoving distance to the redshift z. While 
the intrinsic event rate is still uncertain, we can 


d^'Atobs, (22) 


Figure 8. The same as Figure [3 but the flux is calculated as¬ 
suming that the Ep-E^^iao correlation holds. 


roughly estimate it by using the formation rate of SMSs 
or DGBHs, w hich are theoretically studied in the pre¬ 
vious studies (lAgarwal et al.ll2012t iDiikstra et al.l 12014 
lYue et al.ll20l4). 

lYue et al.l (|2014li studied the formation rate of DGBHs 
in the early Universe. They found that DGBHs can 
be formed from z = 20 to 13 (corresponding to ^ 150 
Myr), and that the comoving mass density of DGBHs 
can be pdcbh 2 x 10® Mq Mpc“®. Since they assumed 
the typical mass of a DGBH as Mdcbh ~ 10® Mq, the 
comoving number density of DGBHs can be evaluated 
as udcbh Pdcbh/-Mdcbh ^ 2Mpc“®[3 Then, we 
can obtain the intrinsic rate of ULGRBs as 'I'grb ~ 
2Mpc“®/150Myr ~ 10“® yr“^ Mpc“®. It should be 
noted that in this rough estimate, we assume that all 
the SMSs collapse to DGBHs after they contribute to 
ULGRBs. This rate may be optimistic. 

Substituting the above value into Equation (|2^ . the 
event rate of the ULGRBs on the whole sky can be ob¬ 
tained as 


AN 

Atobs 


rz=20 


^GKB{z')^T^cr{z'f 


/ 2=13 


dz' 


^6 X 10" 


^GRB _ 

10“8 yr-1 Mpc“® 


yr 


(23) 


The detection rate of the ULGRBs is reduced by the 
beaming factor, Dbeam := 6*^/2 — 3.8x 10“®(6l/5°)^, since 
the off-axis bursts are not detectable. By multiplying 
the beaming factor to Equation (|23)) . the detection rate 
is about one event per year. 

The emission from an expanding cocoon fireball might 
play a key role to raise the detection rate of the event. 
After the jet breakout, the cocoon also breaks out the 
star and evolves like a non-relativistic fireball outsid e 
the star (|Kashivama et al.l l2013t iNakauchi et al.l l2013f) . 
The cocoon emission will be isotropic and free from the 
beaming effect. 


This number density seems lar ger than that of the observed 
typical galaxies. In lYue et al.l II2014I1 . however, they discussed that 
only a fraction of DCBHs can grow up by the mass accretion and 
that most DCBHs do not acquire sufficient mass and escape our 
observation. 
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5.3. Feedback Effects on the Surrounding 
Environments 

Various feedback effects are expected from the super- 
massive collapsars, since they release a huge amount of 
energy. In fact, the total energy of the cocoon fire¬ 
ball discussed above could be as large as Ec ~ 10^®- 
10®® erg, which can be calculated from Equation (TTOll . 
Ec = Ec{t = tb)- Then, the emission from the cocoon 
fireball might be observed as the most energetic super¬ 
nova explosions in the Universe. Moreover, such a vi¬ 
olent explosion could disrupt the host halo, and hinder 
the remnant massive BHs from growing up to SMBHs 
within < 1 Gyr after the BH formation. 

In addition, if heavy elements are produced in the jet 
head and cocoon, they could contribute to the chemical 
enrichment of the host halo. The metal polluted gas will 
induce the formation of second generation of stars. The 
line features in the cocoon emission could also tell us the 
abundance pattern of the nucleosynthesis to confirm the 
SMS origin, although the line may be broad due to the 
high expansi on velocity. 

Rec ently, l.lohnson et al.l (l2013| l and I Whalen et al.l 
(|2013ll considered a very energetic supernova explosion 
of ^ 10®® erg in the first galaxies, and calculated the 
dynamical evolution of the blast wave within the host. 
They found that whereas the blast wave engulfs the en¬ 
tire galaxy, most of its energy is radiated away via effi¬ 
cient cooling processes, so that the swept up matters 
10^ Mq) could recollapse to the host ^ 70Myr after the 
explosion. Thus, such an energetic explosion might not 
hinder the remnant massive BH from becoming super- 
massive within < IGyr after its formation. On the other 
hand, the momentum conservation suggests that the 
host galaxy is expelled if the explosion energy is larger 
than~ lO®®(Mhaio/lO'^M0)(uesc/lOkms“^)(/3c/O.3) erg. 
Thus more detailed calculations are worth while. 

Even after contributing to the prompt emission, the 
relativistic jet has a huge amount of kinetic energy 
^'k.iso = £^ 7 .iso(l - -- 10®'^erg. This can 

leads to bright afterglow emissions at various wave¬ 
lengths (|Toma et al.ll201ll : lloka fc MeszarosI 1200511 . The 
detection of such afterglow emissions could provide us 
rich information about the surrounding environments of 
SMSs such as the density and the chemical composition 
of the circumstellar medium. 

6 . SUMMARY AND GONCLUSIONS 

We investigated whether in the early Universe SMSs 
are able to produce GRBs according to the collapsar sce¬ 
nario. Since SMSs have radii at least as large as RSGs, 
naively it would seem difficult for a relativistic to jet 
reach the surface before the jet engine dies. Actually cal¬ 
culating the jet propagation in SMSs, however, we find 
that jets are able to break out the thick envelope of SMSs. 
This is because the envelope of SMSs is dominated by ra¬ 
diation pressure and has a steeper density gradient than 
RSGs in which the gas pressure dominates. Our conclu¬ 
sion is that SMSs forming in protogalaxies can produce 
violent GRBs. 

Based on empirical rules, we find that the collapse of 
SMSs may be observed as ULGRBs with a duration of 
^ 10^-10® s by the current detectors like BAT. Our op¬ 
timistic estimates indicate rates of detectable GRBs of 


about one event per year. Gomparing with observations, 
we can impose some conditions on the intrinsic event 
rate which is related to the SMS or DCBH formation 
rate. Since GRBs are collimated bursts, beaming re¬ 
duces the rate of detectable events. Therefore, it is im¬ 
portant to consider the isotropic emission accompanying 
GRBs, e.g., the cocoon fireball emissions. Studying the 
detectability of cocoon fireball emission is an interesting 
future work. 

The SMS GRBs are very energetic explosions releas¬ 
ing more than 10®®-10®® erg and sweep up or blow off the 
matter in protogalaxies. As a result, GRBs strongly in¬ 
fluence the mass accretion onto the newborn seed BHs. 
This negative feedback needs to be taken into account 
when studying the growth of remnant BHs. 
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